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Abstract 

We present a non-perturbative determination of the renormahzation factor Zs of the 
scalar density in quenched QCD with overlap fermions. Results are obtained at four 
values of the lattice spacing. By combining with results for the low-energy constant 
E we are able to compute the renormahzation group invariant scalar condensate S in the 
continuum limit with a total accuracy of 7%, excluding dynamical quark effects. Our 
result translates to 5]j^(2 GeV) = (285 it 9MeV)^ if the scale is set by the kaon decay 
constant. We have also performed scaling studies of the pseudoscalar decay constant 
and the vector mass. Our results indicate that quantities computed using overlap quarks 
exhibit excellent scaling behaviour, with small residual lattice artifacts. 
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1 Introduction 



The expectation value of the scalar density at vanishing quark mass, commonly named 
the quark or chiral condensate, plays a central role in QCD at low energies. Spontaneous 
chiral symmetry breaking is signalled by the formation of a non-vanishing condensate, 
and an accurate determination of its value is of great practical interest. Lattice sim- 
ulations of QCD appear well suited for this task, but in order to guarantee a reliable 
error estimation, it is crucial to have control over systematic effects. In particular, to 
ensure that the quark condensate approaches the continuum limit as a power series in 
the lattice spacing a, the renormalization of the bare scalar density must be known with 
good accuracy. It is well known that renormalization factors computed in perturbation 
theory at one loop are not reliable. Further complications arise if the lattice formula- 
tion breaks chiral symmetry explicitly. For instance, in the case of Wilson fermions, a 
cubically divergent term must be subtracted before multiplicative renormalization can 
be applied [1]. 

In this paper we report on a non-perturbative calculation of the renormalization 
factor Zs of the scalar density, using the overlap operator [2] as our fermionic dis- 
cretization in the quenched approximation. We employ the method proposed in [3] and 
compute Zs at four different values of the lattice spacing, ranging from a ~ 0.12 fm to 
0.075 fm. By identifying the bare condensate with the low-energy constant S, which 
appears in effective low-energy descriptions of QCD, we can compute the renormalized 
quantity, given results for S at the corresponding values of the bare coupling in the 
quenched theory. 

Our analysis of the scaling properties of the renormalized condensate indicates the 
presence of only very small cutoff effects of order a? , provided that the non-perturbative 
estimates for the renormalization factor Zs are used throughout. Thus, an extrapolation 
to the continuum limit can be performed in a controlled way. 

Moreover, we have extended the scaling analysis to other quantities, such as the 
pseudoscalar meson decay constant and the mass in the vector channel. In all cases we 
observe an excellent scaling behaviour, with leading cutoff effects of order a^, and thus 
consistent with expectation. To our knowledge, these results represent the first detailed 
scaling study for overlap fermions. 

Results for the quark condensate have already been published by a number of au- 
thors [4-13]. The novelty in this paper is the extension of previous simulations with 
overlap fermions [5, 7, 8] to considerably finer lattice spacings, as well as the strict ap- 
plication of non-perturbative renormalization, enabling us to take the continuum limit. 
Overlap fermions, despite their larger numerical cost, have clear conceptual advantages 
when it comes to studying the problem of chiral symmetry breaking, which is encoded 
in the value for the quark condensate. We stress, though, that our results are valid 
for quenched QCD, and thus great care must be taken if they are to be interpreted in 
the context of the full theory. In particular, the chiral condensate is ill-defined in the 
quenched approximation [14]. 
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2 Renormalization conditions 



Here we briefly recafl the conditions that fix the renormahzation of the scalar and 

pseudoscalar densities in simulations using fermionic discretizations that satisfy the 
Ginsparg- Wilson relation [15,16]. Full details can be found in refs. [3,9]. 

If the regularization preserves chiral symmetry, then the chiral Ward identities 
imply that 

Zs = Zp = XjZ^. (2.1) 

The renormalization factor .Zs, which relates the bare scalar density to the renormal- 
ization group invariant (RGI) density, can then be defined by [3] 



^^^^^^^(rom)(,o) 



(2.2) 

(romp)2=a;i.ef 



In this expression Uyi denotes the RGI quark mass in the continuum limit, in units of 
the hadronic radius ro [17], while m is the bare quark mass that appears in the lattice 
Dirac operator satisfying the Ginsparg- Wilson relation. The expression on the right is 
evaluated at a given reference value, a^rcf , of the square of the pseudoscalar meson mass 
in units of tq. A convenient choice, which we also adopt here, is Xrgf = 1.5736. For 
ro = 0.5 fm this corresponds to mp = mx = 495 MeV. The original data required for 
the determination of IJyi were published in [18], and in eq. (3.1) of [3] IJyi is listed for 
several choices of Xref- 

Since Zg = Zp an alternative renormalization condition can be formulated in terms 
of the matrix element of the pseudoscalar density. If we introduce the shorthand nota- 
tion 

grbare ^ (o]P«(0)lPS), P\x) = (^X'^^'il^){x), (2.3) 



where A" is some flavour matrix, then Zp can be defined via 



Zp = 



(ro^Gr^)(5o) 



(2.4) 



The universal factor Up denotes the RGI matrix element of the pseudoscalar density in 
the continuum limit. Its value can be determined, for instance, using 0(a) improved 
Wilson fermions, and the results presented in refs. [18, 19] then yield 

i7p = 1.802(42) at (ro mp)^ = 1.5736. (2.5) 

In order to compute or Zp, it is clear from eqs. (2.2) and (2.4) that the main task is 
the determination of the value of the bare quark mass, m, and the matrix element G^^^ 
at the point where (ro mp)^ = Xj-ef , for a fermionic discretization based on the overlap 
operator. 
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3 Numerical simulations 



In our simulations we have computed mesonic two-point correlation functions in the 
pseudoscalar and vector channels. We have used the massive overlap operator D^, 
defined by [2] 

Drr, = {l-\am)D + m, D = ^(l--^=^, (3.1) 

where 

A = l + s-aD^, a = -^, \s\<l, (3.2) 

1 ~|~ s 

and is the Wilson-Dirac operator. 

The calculation of the quark propagator proceeds as usual by solving 

Dmi^ = V (3.3) 

for a source field r]. As pointed out in [20], the determination of both chiralities of the 
solution i/j requires some care in the presence of zero modes of the massless operator D, 
especially as the quark mass becomes small. To separate off the zero mode contribution 
we have implemented the strategy outlined in section 7 of ref . [20] , which we briefly 
review here. To this end we shall consider a gauge configuration which has a number of 
zero modes with positive chirality. 

The solution to eq. (3.3) with negative chirality is given by 

P.^P = {DlDm)-^P_Dlr^, (3.4) 

and thus the inversion of I?^Z?„i takes place in the chirality sector that does not contain 
zero modes. The components with positive chirality are obtained from 

= -PoP+r^ + {P+DmP+rH{l - Po)P+v - P-DmP-A, (3.5) 
ml- J 

where Pq is a projector onto the subspace spanned by the zero modes, and whose calcu- 
lation is described in [20]. When implemented in a computer program, eq. (3.5) offers 
complete control over the zero mode contribution. It is also clear that the necessary 
inversion of {P^D„iP^) is performed on a source where all zero mode contributions have 
been projected out. The roles of the positive and negative chirality sectors are obviously 
reversed in the above expressions if the zero modes have negative chirality. 

In our programs we compute P-ijj and P+ijj using the Generalized Minimum Resid- 
ual (GMRES) algorithm [21], which also allows for an inversion of D^n itself. To speed up 
the inversion we have incorporated "low-mode preconditioning" , a technique designed to 
protect against numerical instabilities caused by very small eigenvalues of D\^Dm [20]. 
As we shall see later, the quark masses considered in this work are relatively large and 
hence provide an infrared cutoff, but we found that the inversion can nevertheless be 
accelerated in this way. The presence of zero modes in conjunction with the fact that the 
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low (non-zero) modes arc only known with a certain numerical accuracy requires some 
care in the implementation of low-mode preconditioning for the solution in eq. (3.5). 
Details will be described elsewhere [22]. 

Since the goal of our study is the computation of the renormalized condensate in 
the continuum limit, we have chosen our simulation parameters to coincide with those of 
previous determinations of the bare condensate. To this end we have identified the latter 
with the parameter S computed by matching the spectrum of low-lying eigenvalues of D 
to the predictions of Random Matrix Theory [23] . More precisely, we have concentrated 
on the dataset labeled "B" in that reference, which comprises three different lattice 
spacings at a fixed box size of L = 1.49 fm. We note that a spatial volume of this size 
is sufficiently large to avoid large finite volume effects for masses and decay constants 
at mp mK- In order to improve the accuracy of the continuum extrapolation we 
added a fourth /3-value, (3 = 5.9256, tuned to reproduce the same physical box size for 
L/a = 14. Following the same procedure as in [23], we have determined the low- lying 
spectrum of the Dirac operator and extracted the parameter S. 

The computation of fermionic two-point functions proceeded by setting T = 2L, to 
control the exponential decay of the correlation function in a more reliable way. Our 
simulation parameters are listed in Table 1. As in ref. [23], the parameter s in the 
definition of the overlap operator (c.f. eq. (3.2)) was set to s = 0.4. 



f3 


L/a 


ro/a 


a [fm] 


#cfgs 


5.8458 


12 


4.026 


0.124 


200 


5.9256 


14 


4.697 


0.106 


174 


6.0000 


16 


5.368 


0.093 


200 


6.1366 


20 


6.710 


0.075 


100 



Table 1: Simulation parameters for the determination of mesonic two-point functions 

At each value of the bare coupling we computed quark propagators for three bare 
masses straddling the reference point corresponding to x^f = (romp)^ = 1.5736. We 
added a fourth, heavier value at all but the finest lattice spacing we considered, to study 
the quark mass dependence of mesonic quantities in more detail. Since the quark masses 
here are relatively large, the low-lying spectrum of D cannot induce large fluctuations in 
correlation functions like those observed in the so-called e-regime [24, 25] . Therefore, we 
did not apply the method known as low-mode averaging [25, 26] to enhance the signal. 

In the pseudoscalar channel we used both the left-handed axial current and the 
pseudoscalar density P as interpolating operators, i.e. 

Mx) = (iPr^^P-i^sKx), P{X) = (V^,75V'.)(X), (3.6) 

where P± = ^(1 ± 75), and r, s denote flavour labels. Choosing r ^ s, both of these 
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composite fields were then combined into non-singlet two-point correlation functions 



Cqr(xo) = {Q{x)R{0)) , Q,R= Jo, P. (3.7) 

X 

The correlation function Cjj involves only the left-handed quark propagator such that 
zero modes cannot contribute. By contrast, Cpp includes components of the quark 
propagator whose chirality coincides with that of the zero modes (if any). The latter 
can be separated off by implementing the expression in eq. (3.5). 

The pseudoscalar mass and decay constant, as well as the matrix element Gp^^ 
were extracted from single-cosh fits, after averaging the correlators over the forward 
and backward halves of the lattice. Good plateaus were observed, which served as a 
guideline for choosing our fit intervals. We also computed the current quark mass, 
JTipcAC from 

""^^^^^ = 2 CM^) ' ('-'^ 

where Oq, Oq denote the forward and backward lattice derivatives, respectively. In order 
to compute meson masses in the vector channel, we have considered the two-point 
correlator 

3 

Cvyixo) = a^J2J2i^kix)Vkm, Vk{x) = {^,jM{x), k = 1,2,3. (3.9) 

X k=l 

It turned out to be impossible, however, to obtain a stable plateau for the effective mass, 
by simply using a local source vector in the inversion step. Therefore we applied Jacobi 
smearing on the source 7], as described in [27]. The parameters were chosen such that 
the rms. smearing radius in units of ro was kept constant at approximately 0.6. With 
this choice we were able to improve the stability of the plateau in the vector channel 
considerably. 



4 Determination of renormalization factors 

Our results for masses and matrix elements are summarized in Table 2. 

In order to compute Zg according to eq. (2.2) we have to determine the quark mass 
at the reference point in units of tq. In Fig. 1 we have plotted (romp)^ as a function of 
(rom) at all four /^-values. As can be seen, the data are easily fitted by straight lines, 
but a non-zero intercept is found at all but the largest value of /3: the pseudoscalar 
mass at zero bare quark mass differs from zero by 1 — 2 standard deviations. Since the 
correlation function of the left-handed axial current is free from contributions of zero 
modes, they cannot be responsible for the non-zero intercept. We note however that the 
chiral fits yield x^/dof below 1 even if the extrapolation is forced through the origin. 

By performing local interpolations to the reference point using the three nearest 
data points and subsequently applying eq. (2.2), we obtain the values of Zs, which 
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6.0000 


0.030 


0.01927(7) 


0.217(6) 


0.0346(7) 


0.0814(42) 


0.424(15) 




0.040 


0.02576(7) 


0.247(5) 


0.0356(6) 


0.0612(27) 


0.445(11) 




0.050 


0.03229(7) 


0.273(5) 


0.0366(6) 


0.0501(20) 


0.462(9) 




0.085 


0.05543(9) 


0.352(3) 


0.0403(5) 


0.0342(10) 




6.1366 


0.024 


0.01638(6) 


0.168(5) 


0.0296(7) 


0.0447(21) 


0.360(28) 




0.032 


0.02185(6) 


0.195(4) 


0.0301(6) 


0.0356(14) 


0.378(20) 




0.010 


0.02734(G) 


0.218(4) 


0.0309(6) 


0.0305(11) 


0.389(15) 



Table 2: Results for meson masses and decay constants computed at several values of quark masses at 
each lattice spacing. The results for amp and aFp'""'* were extracted from correlators of the left-handed 
a:xial current. 

are tabulated at each /3-value in Table 3. The typical accuracy of our determination 

is around 5%. It should be noted that the precision is partly limited by the accuracy 
of the published value of Um, which is about 3% [18]. We estimate that pushing the 
precision of our determination of Zs to that level would require a four-fold increase in 
statistics. 

In Fig. 2 we plot our results for Zs versus /3. It has become customary to rep- 
resent results for renormalization factors at different values of the bare coupling by 
interpolating curves. Using a simple polynomial ansatz in (/3 — 6) yields 

Zs(/3) = 1.045-0.899(/3-6)+4.36(/3-6)^ s = 0.4. (4.1) 

This formula describes Zs with an estimated error of 5% in the studied range of P, i.e. 
5.8458 < P < 6.1366. We emphasize that our determination is valid only for the case 
s = 0.4 in the definition of the Neuberger-Dirac operator, eqs. (3.1) and (3.2). 
The perturbative expression for Zg at one loop is 

zrim) = [i + g'o ln(a/.) + + Oia'o)] , (4.2) 

where = 0.147107 for our choice of s = 0.4 [28, 29]. The factor rnj^{ii)/M was 
computed previously in [18]. The mean-field improved version of reads [3] 



6 



(roiTLp) 

3.5 
3 

2.5 
2 

1.5 
1 

0.5 




/S = 6.1366 
• /S = 6.0000 

/S = 5.9256 
. /S = 5.8458 




0.1 0.2 0.3 0.4 



Figure 1: (romp) as a function of rom. The horizontal dashed line represents the reference point 
(romp)2 = 1.5736. 









Za 


5.8458 


1.28(6) 


1.33(4) 


1.710(5) 


5.9256 


1.19(7) 


1.20(4) 


1.611(3) 


6.0000 


1.05(5) 


0.88(6) 


1.553(2) 


6.1366 


1.01(4) 


1.02(5) 


1.478(2) 



Table 3: Non-perturbative determinations of Zs, Zp and Za- 



Zfigo 



m- 



ifi) fl + s 



M 



l + s 



ln(a/i) + + 



'■0 



l + s 



(4.3) 



where = Qq/uq is the boosted coupling, s = 3 + {s — 3)/uq, with Uq being the average 
plaquette. The comparison of our numerical results for Zs with perturbation theory is 
shown in Fig. 2. The mean-field improved perturbative expansion comes quite close to 
the non-perturbatively determined values for /? ^ 6.0 but falls short by more than 20% 
below P = 6.0. Unsurprisingly, perturbation theory in the bare coupling fares a lot 
worse in the entire range of couplings studied here. 

The results for Zp, computed according to eq. (2.4), are listed alongside those for 
Zs in Table 3. The renormalization conditions for Zo, and Zp imply that the two must 
be identical up to effects of order a^. Indeed, we observe hardly any difference at our 
level of accuracy, except at /9 = 6.0. In our view, the most likely explanation for this 
deviation is a statistical fluctuation. 

In order to include the pseudoscalar decay constant in the scaling tests described 
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Figure 2: Zs as a function of /3. The solid line 
denotes the fit of eq. (4.1). The dotted and dashed 
curves represent the results of bare and mean-field 
improved perturbation theory at one loop order. 



Figure 3: The quark mass dependence of — — — . 

o ^ J' mpcAC 

The value of [3 increases from top to bottom. Za 
is defined as the value of this ratio in the limit of 
vanishing quark mass. 



below we also computed the renormalization factor of the axial current, Za- Using the 
PCAC relation and Z^a = l/-^p one can define 

77? 

Za = lim . (4.4) 

m^O mpcAC 

We found the ratio m/nipcAC to depend only weakly on the bare mass (c.f. Fig. 3). Za 
could then be determined by extrapolating m/mpcAC linearly in m to the chiral limit. 

5 The renormalized condensate 

Having determined the renormalization factor of the scalar density in a range of bare 
couplings, we can now compute the renormalized condensate in the continuum limit, by 
combining the results for Z^ with estimates of the bare condensate. 

In effective low-energy descriptions of QCD with A'^f = 3 quark flavours, the quark 
condensate is identified with the low-energy constant S via 

- (V^V-) = S. (5.1) 

In the quenched theory, however, the condensate —{^l^ij)) is not defined, owing to the 
presence of infrared divergencies as the chiral limit is approached [14]. Nevertheless, the 
low-energy constant S can be determined in quenched QCD, for instance, by comparing 
lattice data of suitable quantities to expressions of Chiral Perturbation Theory or chiral 
Random Matrix Theory. Although in this case the identification of S with the quark 
condensate is rather dubious, we shall nevertheless proceed to compute a renormalized 
"condensate" , by assuming that estimates of S in the quenched theory renormalize like 
the scalar density. 
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Our input quantities are thus the renormahzation factors Zs of Table 3 and results 
for E, determined by matching the low- lying eigenvalues of the Dirac operator in the 
e-regime to the predictions of the chiral unitary random matrix model according to [23] 



k = l,2, 



(5.2) 



Here, {Xk)u is the expectation value of the fcth eigenvalue in the topological sector 
with index u, and Cfc denotes the kth scaled eigenvalue in the matrix model. In ref. 
[23] it was found that good agreement with random matrix behaviour is observed for 
lattice volumes V of at least (1.5 fm)^. In other words, the value of S extracted from 
eq. (5.2) depends neither on the particular eigenvalue, nor on the topological sector, 
within statistical errors. 

Using the results for S from Table 3 of [23] (i.e. the runs labelled Bo,Bi and B2), 
supplemented by our data at /3 = 5.9256, we plot the renormahzation group invariant 
condensate S in units of tq versus (a/rg)^ in Fig. 4. If the non-perturbative estimates 
for Zs are used, the results for TqE show a remarkably flat behaviour, which not only 
indicates small residual cutoff effects, but is also consistent with the expectation that 
the leading lattice artefacts of our fermionic discretization should be of order a^. 

Figure 4 also reveals that employing mean-field improved perturbation theory for 
Zs produces a significant slope in rgS as the continuum limit is approached. Although 
this procedure apparently yields a consistent value of TqE in the continuum limit, it is 
equally obvious that the perturbatively renormalized result serves as a poor estimate 
for the condensate at non-zero lattice spacing. 
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Figure 4: Continuum extrapolation of TqE. Full 
circles denote the results obtained using non- 
perturbative renormalization factors, while open 
squares represent values resulting from applying 
mean-field improved perturbation theory. 



Figure 5: The variation of TqE in the continuum 
limit, arising from choosing different eigenvalues 
and topological sectors in the determination of the 
bare condensate. The solid and dashed lines rep- 
resent the result for k = 2, = 1 which is used 
for our main result. 



Our results for TqE at all values of (3 and in the continuum limit are listed in 
Table 4. Here we have used E as determined from {Xk)u for k = 2 and \u\ = 1. We 
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note that the variation in the value of r^S from choosing different A^'s and topological 
sectors is well within the statistical fluctuations after taking the continuum limit. This 
is illustrated in Fig. 5, where we plot the continuum results for all possible choices of 
Afe and We emphasize that this variation should not be regarded as a systematic 
uncertainty, since all choices are equivalent, if random matrix theory does indeed give 
an accurate description of the low-lying eigenvalues, and hence we refrain from quoting 
an additional error. 
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0.293(21) 


0.294(9) 


2.32(29) 



Table 4: Renormalization group invariant quark condensate, kaon decay constant and K^-mass, in 
units of ro- 

Our result in the continuum limit is thus 

r^E = 0.293 ± 0.021 (5.3) 

for the renormalization group invariant condensate. In the MS-scheme at 2GeV we 
obtain after division by mj;jg(2 GeV)/M = 0.72076 [18] the value 

rgS|^(2 GeV) = 0.406 ± 0.029. (5.4) 

These are the main results of our calculation. To our knowledge, these are the first 
estimates of a quantity in the continuum limit, computed using overlap fermions. We 
emphasize that the quoted errors include all uncertainties, except those due to quench- 
ing. 

As is well known, the calibration of the lattice spacing is ambiguous in the quenched 
approximation, and thus any conversion into physical units is only illustrative. Here we 
perform such a conversion using either the kaon decay constant or the nucleon mass to 
set the scale. Ref. [18] quotes 

roFK\/2 = 0.415 ± 0.009, Fk = 113 MeV, (5.5) 

in the continuum limit, while a continuum extrapolation of the nucleon mass data of [30] 
in units of ro yields 

romN = 2.670 ± 0.042, rriN = 939.6 MeV. (5.6) 
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For the condensate in the MS-scheme at 2 GeV we then obtain 

v_r9P - / (2^^ ± 9MeV)3, scale set by Fk , 
^Msi^'-evj - j ^26l±8MeV)3, scale set by mN ' ^ ' 

These findings are consistent with previous observations that the typical scale ambi- 
guity for a quantity with mass dimension equal to one is of the order of 10%. Recent 
calculations of the renormalized condensate [4-13] yield similar values compared to our 
results. 



6 Further scaling tests 

The leading cutoff effects of fermionic discretizations based on the Ginsparg- Wilson 
relation are expected be of order a^, and indeed, this expectation has been confirmed in 
our scaling study of the quark condensate. In this section we shall extend our analysis 
of cutoff effects to quantities like the pseudoscalar decay constant and the meson mass 
in the vector channel. 




0.01 0.02 0.03 0.04 0.05 0.06 




0.01 0.02 0.03 0.04 0.05 0.06 



(a/rj= 



Figure 6: Continuum extrapolation of roJk. 
Full circles denote our results, while the open 
squares are the data of [18], employing 0(a) im- 
proved Wilson fermions. The full triangles are our 
data with Za from mean-field improved perturba- 
tion theory. 



Figure 7: Scaling behaviour of rorriK* • The 
meaning of the full circles and open squares is as 
in Fig. 6. The open circle results from an alterna- 
tive fit with a fit range of xo/a G [5, 11] instead of 
xo/a e [8, 11] (fuU circle). 



To this end we have assumed that aFp^'^'^ and amy depend linearly on (amp)^ and 
performed a linear interpolation to the point where (romp)^ = (rom-K)^ = 1.5736. Thus, 
our aim is to investigate the scaling behaviour of and mx* . The renormalized kaon 
decay constant is obtained after multiplication with the factor Za listed in Table 3. In 
Table 4 we have compiled the results for roFx and romx* at the various values of (3, as 
well as in the continuum limit. The corresponding continuum extrapolations are plotted 
in Figures 6 and 7. 
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For the kaon decay constant we observe a flat approach to the continuum hmit, 
consistent with a hnear fit in (a/ro)^, provided that the non-perturbative estimate for 
is used. The perturbatively renormaUzed Fk is subject to larger lattice artefacts, 
and the resulting continuum value is roughly consistent. In Fig. 6 we also show the 
continuum extrapolation of the same quantity from ref. [18], where tq-Fk was computed 
using 0(a) improved Wilson fermions. In the continuum limit our data agree remarkably 
well with those of ref. [18], but for overlap fermions the residual cutoff effects at lattice 
spacings of around 0.1 fm, i.e. at (a/ro)^ ~ 0.035, are apparently much smaller. 

The scaling behaviour of the K* mass is also flat, except at our coarsest lattice 
spacing. A closer inspection of our fits to the two-point function shows that the value of 
amy at P = 5.8458 depends strongly on the chosen fit range. Extending the fit interval 
to smaller timeslices leads to a significant increase in the value of rorriK* , as indicated 
in Fig. 7. Owing to the uncertainty in the value of roiriK* as a result of using different 
fit intervals, we exclude the coarsest lattice from the continuum extrapolation, despite 
the fact that the alternative result is apparently consistent with a linear behaviour up 
to (a/ro) ~ 0.06. Nevertheless we also confirm good scaling behaviour for the vector 
mass; as our values for /? > 5.8458 are mutually consistent with each other, as well as 
with the results of ref. [18]. 

7 Conclusions 

We have presented the first comprehensive scaling study of quantities computed using 
overlap fermions. A major part of our calculation was devoted to the determination of 
the renormalization factor Zo, of the scalar density. Thereby we were able to present 
a conceptually clean determination of the renormalized low-energy constant S in the 
continuum limit of quenched QCD, with a total accuracy of 7%. 

Besides studying the continuum extrapolation of rgS we also performed scaling 
studies of the pseudoscalar decay constant and the mass in the vector channel. For 
all three quantities computed using overlap quarks we observed an excellent scaling 
behaviour, resulting in a fiat approach to the continuum limit. This is signified by the 
fact that the results in Table 4 at any finite value of P and in the continuum limit are 
practically the same, at least at our level of accuracy. We note, however, that a flat 
continuum behaviour is only observed for S and Fk, if non-perturbative estimates of 
the respective renormalization factors are employed. Our values for roFk and rQiriK* in 
the continuum limit are in very good agreement with those of refs. [18,31]. 

Owing to their good scaling properties, overlap fermions are an attractive discretiza- 
tion for the computation of phenomenologically interesting quantities, despite the large 
numerical effort involved in their simulation. 
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